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Abstract. 

General recurrence relations for arbitrary non-diagonal, radial hydrogenic ma- 
trix elements are derived in Dirac relativistic quantum mechanics. Our approach 
is based on a generalization of the second hypervirial method previously employed 
in the non-relativistic Schrodinger case. A relativistic version of the Pasternack- 
Sternheimer relation is thence obtained in the diagonal {i.e. total angular momentum 
and parity the same) case, from such relation an expression for the relativistic virial 
theorem is deduced. To contribute to the utility of the relations, explicit expressions 
for the radial matrix elements of functions of the form and /3r^ — where /3 is a 
Dirac matrix — are presented. 
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I. Introduction 



Recurrence relations for matrix elements are very useful tools in quantum calcu- 
lations [1-5] since the direct computation of such elements is generally very cumber- 
some. An interesting example is the Blanchard's relation which is a useful recurrence 
formula for arbitrary (i.e. not necessarily diagonal) non-relativistic matrix elements 
of the form {nili\r^\n2l2) , where the \nl) stand for non-relativistic hydrogenic radial 
energy eigenstates; according to this relation, once any three successive matrix ele- 
ments of powers of the radial coordinate, r, are known, any other can be deduced in 
terms of the previously known ones. The Blanchard recurrence relation was derived 
more than twenty five years ago using a calculation-intensive method [6], which is 
not surprising since in general recurrence relations are rather diflicult to obtain. 

Trying to overcome such difficulties, different approaches have been proposed for 
obtaining general recurrence relations; some of them are based on algebraic methods, 
others use sum rules and hypervirial theorems [4,7-10]. In particular, a hypervirial 
result has been employed to obtain the Blanchard relation in a more compact way 
than in the original deduction [11]. In relativistic quantum mechanics, on the other 
hand, and despite its physical and its possible chemical interest [7,8,12,13], excepting 
for the results reported in [14], there are not as yet general recurrence relations which 
could be used for calculating matrix elements of powers of the radial coordinate in 
terms of previously known ones. We should mention though that there are previous 
efforts in such direction, in which closed forms for certain relativistic matrix elements 
have been evaluated [16,17] and certain, mostly diagonal, recurrence relations for 
relativistic and quasirelativistic states have been calculated [3,9,12,17] . 

In this paper we employ a relativistic calculation inspired on the hypervirial 
method [11] to deduce a recurrence relation for the, in general, non-diagonal radial 
matrix elements of succesive powers of and of /i^r^ — where /3 is a 4x4 Dirac matrix 
[18-20] — for relativistic hydrogenic states in the energy basis. The assumptions we 
use are that the nucleus is point-like and fixed in space, and that a description 
using the Dirac equation is valid. We first study the recurrence relations in the 
general case, in which the matrix elements are taken between states with different 
principal quantum numbers rii ^ 77.2,, different total angular momentum quantum 
numbers ji 7^ j2, mj^ 7^ mj2, and, as we use the quantum number e = (— 
instead of parity for labelling the hydrogenic eigenstates, where ei ^ 62- We find 
that in general the recurrence relations depend on matrix elements of both powers 
of r and of /3r. In practical terms this means that we need two recurrence relations 
as the relativistic version of the single-equation Blanchard relation [Eqs. (6) and (7) 
of section II]. Given its special interest, we study in particular the case where the 
total angular momentum and parity become equal, ji = j2 and ei = 62, in the two 
states — not mattering the relative values of the principal quantum number n. We 
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also address the completely diagonal case where ni = n2, ji = 32-, and e\ = €2- 
Both of the particular cases mentioned above require special treatment for avoiding 
possible divisions by zero in the general expressions; such results are immediately 
used to obtain a relativistic version of the Pasternack-Sternheimer rule [21] and to 
obtain an expression for the relativistic virial theorem [9,14,22]. 

This paper is organized as follows. In section II we review the second hypervirial 
scheme used for deriving the non-relativistic Blanchard relation. In section III, after 
obtaining the radial Hamiltonian useful for implementing the hypervirial result in 
relativistic quantum mechanics, we proceed to use it to deduce by a long but direct 
calculation the relativistic recurrence formulae. In section IV we study in particular 
the diagonal case (72 = ji, £2 = ei) to derive the relativistic Pasternak-Sternheimer 
rule and use it (when ni = n2) to obtain a version of the relativistic virial theorem. 
In the Appendix, we obtain explicit expressions for diagonal and non-diagonal matrix 
elements for any power of r and of I3r between radial relativistic hydrogenic states. 
As it becomes evident, such results are rather cumbersome for relatively large values 
of the power; for small values, on the other hand, they are better regarded as starting 
values for the recurrence relations derived in section III of this article. Furthermore, 
these results can be of utility not only for relativistic atomic and molecular studies 
but also for evaluating matrix elements of interactions designed to test Lorentz and 
CPT invariance in hydrogen [14,29]. 

II. The non-relativistic recurrence relation 

Both the Blanchard relation and its predecesor the Kramers selection rule, were 
originally obtained employing directly the Schrodinger equation together with ap- 
propriate boundary conditions and, at least in the former case, a great deal of com- 
putations [1,6]. A much simpler approach based on a generalized hypervirial result 
and certain Hamiltonian identities has been developed to simplify the computations 
leading to the Blanchard relation [11]. This technique seemed to us an appropriate 
starting point for deriving relativistic recurrence formulae. It is with such relativis- 
tic extension in mind that we review in this section the hypervirial method as it 
is applied in non-relativistic quantum mechanics. In this section, as in most of the 
paper, we employ atomic units h = m = e = 1. 

The idea is to start with the radial Schrodinger equation for a central potential 
V{r) written in the form 

Hk\nklk) = En^i^\nklk), (1) 

where \nklk) = i^nkh (^) h energy eigenfunction and its corresponding 

energy eigenvalue with principal and angular momentum quantum numbers, Uk and 
Ik, respectively; k is just a label, and Hk-, the non-relativistic radial Hamiltonian, is 
given by 
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Although we want to calculate the radial matrix elements of terms of the form 
r^, it is best for our purposes to consider first matrix elements of an arbitrary radial 
function /(r). With such choice we can readily show [11] that 



{Ei - Ek){ni k \f{r)\nk k) = k \ [-\f" " / ^ " + k), (3) 

where we use A^^ = /j (li + 1) — k {Ik + 1), = E^^i^, and the primes stand 
for radial derivatives. Please recall that the matrix element of an arbitrary radial 
function /(r) is 

/"OO 

{nik\f{r)\nkh) = / rV^, JO /(r)Vn.;. (r)dr. (4) 

To establish the result we are after, we apply the previous result (3) to the radial 
function ^(r) = Hif{r) — f{r)Hk, to find 



2{Ei-Ekf{nih\f{r)\nkh) = 

{m k I {Hi {HJ{r) - f{r)Hk) - (HJir) - f{r)Hu)Hu+ (5) 
Hi (HJir) - f{r)Hk) - {HJ{r) - f{r)Hk)Hu) W h)- 

This is the generalized second hypervirial valid for arbitrary radial potential energy 
functions, V{r), introduced in Eq. (8) of Ref. 11. 

The second hypervirial takes a particularly simple form when /(r) is a power 
of the position, let us say /(r) = r^"*"^; using this expression for /(r) and restricting 
ourselves to the Coulomb potential, V{r) = —Z/r, we obtain [11], after a long — but 
much shorter than in [6] — direct calculation, the Blanchard relation 



X{Ei-Ek) {nili\r^'^'^\nklk) = co{nili\r^\nklk) + ci{nili\r^ ^Irikh) 

I r I A-2, , N 

+ C2{nili\r \nklk); 
where the hydrogenic energy eigenvalues are Ea = —Z^/2n'^, independent of and 



{k-lk){k + lk + l) 



nt 



+ A(A + 2) ^ + 



nt 



ci = -2ZA(A + 2)(2A + 1) 



(7) 



C2 = - (A + 2) [A^ - {Ik - kf] [{Ik + + 1)' - A^] . 
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From this result we can also obtain, as special cases of the Blanchard recurrence 
relation (6), first the Pasternack-Sternheimer selection rule [21]: 

{nih l^lrikk) = 0, (8) 

saying that the matrix element of the potential vanishes between radial states of 
central potentials when their angular momenta coincide and when the corresponding 
energy eigenvalues depend on the principal quantum number only. Second, in the 
completely diagonal case {i.e. rii = n/c, li = Ik), we can further obtain the non- 
relativistic quantum virial theorem [9] 

{V) = -Z{^)^2{E). (9) 
r 

As we exhibit in section IV, we can obtain analogous results using our recurrence 
relations in relativistic quantum mechanics. 

III. The relativistic recurrence relations. 

In this section we apply the method sketched in section II to the relativistic 
Dirac case. We clearly need to start with a radial Dirac Hamiltonian analogous to 
(2). To obtain such Hamiltonian we start with the Dirac Hamiltonian Hd and the 
corresponding time-independent Dirac equation for a central potential 

HD^ca-p + pc^ + V{r), i?z)*(r) = ^*(r); (10) 

where we are using atomic units, a and P are the 4x4 Dirac matrices [18-20], which 
in the Dirac representation are given by 

o). -O' 

where the I's and O's stand respectively, for 2x2 unit and zero matrices and the 
a is the vector composed by the three 2x2 Pauli matrices cr = (ax, cry,az). Please 
notice that, despite the selection of natural units we shall, where it aids interpreta- 
tion, reinsert the appropriate dimensional factors in certain equations. The energy 
eigenvalues are given explicitly in Eq. (63) of section V. The Hamiltonian Hd is 
rotationally invariant, hence the solutions of the Dirac equation (10) can be written 
in the alternative but entirely equivalent forms [19,23] 



*(r,^,0) = - =- I, (12) 

where x^m^ and x-nm^, or 3^j„j and 3^j^, are spinor spherical harmonics of opposite 
parity, and k = — e(j -|- 1/2) is the eigenvalue of the operator A = (3{1 + S • L) which 
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commutes with Hd (where S = o" ® / = diag(o", a)). The second form in (12) is the 
preferred in Ref. 18. Parity is a good quantum number in the problem because central 
potentials are invariant under reflections; parity varies as (—1)^ and, according to 
the triangle's rule of addition of momenta, the orbital angular momentum is given 
by / = J ± 1/2. But, instead of working directly with parity or with k, we prefer the 
quantum numbers j and e, introduced above, which can be shown also to be 

fl If/=j + |, 

(13) 

[-1 Ifl^j-l, 

thus / = J + e/2 in all cases. We also define = j — e/2; accordingly, the spherical 
spinor depends on / whereas the spherical spinor y'j^, which has the opposite 
parity, depends on I' . Writing the solutions in the form (12) completely solves the 
angular part of the problem. 

To construct the radial Hamiltonian, we use the relation 

(« ■ r)(« ■ p) = (E ■ r)(E ■ p) = r • p + iE • L; (14) 

we then use = [L + (1/2)S]^ = L^ + S- L + 3/4 but for expressing the term 
L • S, we also need an expression for acting on the eigenfunctions (12). Directly 
from this equation we see that when is applied to any central potential state, the 
big component of the state function behaves with the orbital quantum number / = 
j + e/2, whereas the small one does so with the orbital quantum number I' = j — e/2; 
we have then, 

l{l + l)=j{j + l) + e{j + ^) + ^, (15) 
for the big component, and 

ni' + l)=jU + ^)-e{j + l) + l, (16) 

for the small one. The action of upon a solution of the form (12) is therefore 
always of the form 

L'=jU + l) + pe{j + ^) + \, (17) 

where /3 is the Dirac matrix (11). From this result we obtain the term L • E and, 
substituting it into (a • p), we finally obtain 

{a-p) = ar\pr-i/3-{j + l)i (18) 

where 
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ar = -oe ■ r, pr 

r 



dr 



(19) 



We are now ready to write the relativistic radial Hamiltonian, and the corre- 
sponding radial Dirac equation, as 



Hk = car Pr - iP— I Jfc . o 
r \ 2 

Hkipkir) = Ekipkir), 
where we introduced the purely radial eigenfunctions 



(3c^ + V{r), 



in a 2 X 2 representation where, (3 = diag(+l, —1), = 
Dirac equation becomes then [14,19] 



-1 
-1 



(20) 



(21) 

, and the radial 



c'' + (Vkir) - Ek) -ci-ekijk + 1/2) /r + d/dr) 

c{ek [jk + 1/2) /r + d/dr) + {Vk{r) - Ek) 



(r) 



rikJkf- 



kir) 



= 0. (22) 



Though this explicit representation can be used for our problem [24,25], it is not 
really necessary since all our results are representation independent. 

The relativistic recurrence relation we are after, can be deduced using a similar 
reasoning as the used in section II for the non-relativistic case. Let us first calculate 
the non-diagonal matrix element of an arbitrary radial function /(r) 



iE2 - Ei){n2j2e2\f{r)\nijiei) = 

{n2j2e2\H2f{r) - f{r)Hi\ni jid) ^^3) 

= -ic{n2 j2e2\ar (^f'(.r) + ^/5/(^')^ l^u'iei), 

where from now on the labelling in the kets stand for the three quantum numbers 
rifc, jfc, and Ck, we have defined = £2(2^2 + 1) — ei(2ji + 1), and the matrix 
elements of radial functions are calculated as 



{n2j2 €2\f{r)\niji ei) 
{n2j2 e2\Pf{r)\niji ei) 



= / fir) {F^{r)F,{r) + G*,{r)G,{r))dr, 
= J fir) {F*{r)F^{r)-Gl{r)G^{r))dr. 
7 



(24) 



where the subscripts stand for the 3 quantum numbers specifying the state. 

We next proceed to calculate a "second order iteration" by substituting /(r) 
^(r) = H2f{r) — f{r)Hi in the last expression. Let us calculate first and 



^,6^(2^+1)^ f /'(r) + ^Pf{r)] - icar (fir) + ^/3/(r) ) (F(r) 



2r 



2r 



(25) 



and 



^^1 = 4 (/'W - ^mr)) - [fir] - ^mr)) |:+ 

- (/'(O - ^/5/(r)) - ica. (/'(r) + ^/3/(r)) (V^(r) + Pc') . 

(26) 

Then, we write down the difference of the matrix elements associated with Eqs. (25) 
and (26) 

{E2-Eif{n2 32 e2|/(r)|niji ci) = 



c'^f3nr)+c' /(r) + 2zcV/5 (/'(r) + ^/^/w) |niji6i). 

where we have defined A^i = e2(2j2 + 1) + ei(2ji + 1). Please notice that here and 
in what follows we are always assuming A^^ ^ 0. 

This last expression (27) is the direct relativistic equivalent of the generalized 
second hypervirial [Cf. Eq. (5) above]. The expression involves the operator d/dr, 
but here, due to the presence of Dirac matrices in the result, we cannot use the trick 
employed in the non relativistic case where we took advantage of the Hamiltonian 
to simplify the calculation [11]. Instead, let us calculate the following second order 
iteration for non-diagonal matrix elements 



{n2j2e2\H2^ + ^Hi\niji a) = {El - Ef){n2j2 e2\f{r)\ni ji d) = 



4^2 



fir) - 2icar[f'ir) + ^/3/(r)] V(r)^ jiei); 



due to the presence of Dirac matrices in our results, we also require to calculate 
non-diagonal matrix elements for expressions involving arf{r) and (3f{r), namely 



H2 i-iarf{r)) = 
fir) 



nr)-f{r)- 



d , e2 



^(2j2 + l)/3/(r) 



and 



{-iarf{r))H^ = 



cf{r) 



adding up these two last expressions, we get 



ic^ar(3fir) - iarV{r)f{r), 



(29) 



ic^ar^fir) - iarV{r)f{r); 

(30) 



{E2 + Ei){n2j2 €2! - ictrfir)\niji ei) = 

9rf(r') d A~ (31) 

{n2j2 621 - - cf'ir) - 2c/(r)-^ + c-^/3/(r) - 2iarV{r)fir)\n, h ^i)- 

r dr 2r 

From the matrix element of H2 {—iarcf{r)) — {—iarcf{r)) Hi, we can obtain 
{E2 - -Ei)(n2 j2e2| - iarf{r)\niji ei) = 

A+ (32) 
(^2 J2 e2| - cf{r) + c-^/3/(r) + 2cHar^f{r)\mji ei); 

2r 

proceeding in a similar way for H2 {(3f{r)) + {(3f{r)) Hi, we get 



{E2 + Ei){n2 j2 e2\(3f{r)\ni ji ei) =(712^2 e2|ic/3ar/'(r) - icar^-fir) ^^^^ 

+ 2p+/3F(r)] /(r)|niji ei). 

Equations (23-33) are the basic equations of our problem. To proceed, we 

consider, as in the non-relativistic case, radial functions of the form /(r) = and 
insert the explicit expression for the Coulomb potential: V{r) = —Z/r. Let us 
mention though that our results can be generalized to other power of potentials, 
such as the Lennard- Jones potentials [26]. 
Substituting /(r) = in (28), it follows 



[El- El) {n2 32 e2|r^|ni jiei) = 

{n2ke2\c^ [^^-A(A + 1) 



A-2 , 2^^ A-2 _ 2c^\r^-^— + 



2icar ( A + ) r^-^V{r)\ni h ei); 



hence, we can eliminate the term containing the derivative operator in this last 
equation, using /(r) = r^~^ in Eq. (31), to get the result 



{E^ - E^) (712 j2 e2\r^\ni ji d) = 



{n2j2e2\c'^^^r'-' + c'^P{l - X)r^-' - ica^/^A^lr^- V(r)+ (35) 
+ {E2 + Ei)X i-icar) r^- Vi Ji ^i); 

in this last equation, we use Eq. (23) to eliminate the term with —icar 
get 



{E2 - El){n2j2 ^2\r^\niji ei) = (n2j2e2|c^ 



A2lAjl 



(1 - A)/3 



4 2 

2Z [icarr^-'^{l - A) - {E2 - Ei)r^-'^] - {E2 + ^^OAica^r^" Vi ii ei 
Now, from Eq. (32) with /(r) = r'^~^ we get 



,A-2 



(36) 



(^2 - ^i)(n2i2e2| -iarr^ Viii^i) = 



n2j2 £2! - c(A - l)r^"2 + c-^/3r-^"^ + 2ic^Q;r/?mr^~-^|ni ji ei) 



(37) 



and, using /(r) = r in Eq. (33) to ehminate the term 2icar(3mr from the above 
equation, we obtain 



{E2 - Ei){n2j2e2\ - ioCrT^ Vl Jl ^l) = 



{n2 j2 £2! - c (A - 1) + c^pr 



,A-2 



4c3 . 4cZ ^ XI, 
A A 



which can be written as 



J {E2 + E^) Pr^ + ^ (-ia,) A^,r^-'+ 



(38) 



(£^2 ~ -E'l 
4c3 



A" r2 
A 



A 2 



('^2 ^2 e2|(-iarr )|rai ji ei) = (n2 ^'2 £2! - c(A - l)r + 

f|£/3r^-i _ ^(^2 +£;i)/3r^|nu-i ei). 
A A 



(39) 



We can also obtain a new relationship for the matrix elements of —iarr^~^, 
using Eq. (23) with /(r) = r^, and substitute the result in Eq. (37) to eliminate the 
term 2iarPmr^~^ 
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(E2 - Ei){n2j2e2\ - ioirr \nijiei) = 
{n2 32 e2 1 - c (A - 1) r^'^ + c^(5r>^-^ + ^ 



-lar) r 



A-l 



'21 



4c 



-(^2 -£;i)r^|ni ji ei). 



'21 



Rearranging terms, we obtain 



{E2 — El) 



4c2A 



^21 



{n2 j2 e2\{-iarr^ ^)\nijiei) = 



{n2 32 €2! - c(A - l)r 



(40) 



(41) 



The relation we are looking for follows from this last result and Eq. (36). We use 
succesively r^^^ and r'^^^ from Eq. (41) to eliminate the terms 2{E2-\- Ei)Xicarr'^~^ 
and 2zco;r-r'^~^(l — A) that appear in Eq. (36) to finally get [14] 



CQ{n2 32 (^2\r^\ni ii ti) = ^ Ci(n2 J2 e2|r-^ ji ci) + ^ cii(n2 J2 e2|/3r^ *|nijiei), 

(42) 



1=1 



i=2 



where the numbers Cj, i = 0, . . . 3 are given by 

{ El-El){E2-Ei)^^^ 
{E2- Ei)A^^-4c^X ' 
2Z(E2-Ei)2AJi 



CO = 



Cl = 



(S2-£^i)A^i-4c2(A-l) 



C2 = C 



C3 = 



2A2-1A+ _^2^^^_-^^ iEi+E2)A 



21 



4 ' ''(E2-Ei)A2i-4c2A- 
2Zc2(A- 1)(A- 2)A2 



'21 



(£;2-£;i)A2i-4c2(A-i)' 

and the numbers di, i = 2 and 3, by 



d2 = c'^ 



(1-A) + 



A(£;2 + £^i)A 



+ 
21 



{E2 - Ei)A^^ - 4c2A, 



(43) 



(44) 



(E2-£^i)A2i-4c2(A-l) 

As we may have expected, we need to know six matrix coefficients instead of 
only three as in the non-relativistic case. This is a consequence of the fact that, in 
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the Dirac case, we have to deal with the big and the small components in the state 
function, doubling in this sense the "degrees of freedom" of the system. 

It does not seems to be possible to avoid the /9-dependency in Eq. (44), and thus, 
taken on its own, Eq. (41) does not allow the computation of (^2 j2 ^2\^^\^i ji Ci) in 
terms of the {'n2 j2 G2\f^'^~"'\'n'i ji ei) , a = 1,2,3. The situation is not hopeless though 
because it is still possible to obtain another recurrence relation for non-diagonal 
matrix elements of /3r^ simply by eliminating the term — ia^-r^"^ between Eqs. (39) 
and (41). In such a way we get 



eo('^2 j2e2|/3r^|ni ji ei) = 6o(?i2 J2 e2|r^|ni ji ei) + 62(^2 J2 e2|r^ ^|?^iiiei) ^^^^ 
+ ei(n2 j2 e2|/3r^"Viii ^i) + 62(^2^2 e2\(3r^~'^\ni ji ei), 

where the numbers bi and ei i = 1,2,3 are given by 

60 =4A [(^2 -Elf- 4c^] , 
62 =c^(l-A) - 4X'] , 

eo =2(^2 + Ei)[iE2 - Ei)A^i - 4(^X], ^^g^ 
ei =4Z[4c^X-{E2-Ei)A^^], 

e2 =^^[{A-,r-AX% 

Equations (42) and (45) together are the useful recurrence relations in the relativistic 
Dirac case. 

IV. The diagonal case = (j2 = ji, £2 = ei)- 

In the results of the last section we always assume A21 7^ 0, but in order to study 
the diagonal case we must have ei = €2 and ji = J2; this in turn imply A21 = and 
(as always!) = A+ 7^ 0. To deal with the diagonal case we start all over again. 
The equation set for this case is particularly simple, first from Eq. (23) we have 

{E2 - El ) {n2 je\f{r)\nije) = (n2 j e \ {-icarf (r) Im je), (47) 

then we can proceed to calculate the second order iteration by substituting, as in 
the previous section, /(r) ^ ^_ = H2f{r) — f{r)Hi in (47) to obtain 



iE2-E,f{n2je\f{r)\mje) = {n2je\-c'nr) + c'^f'{r)^+ ^^^^ 

2ic^arPf'ir)\nije); 

and then substitute /(r) — > ^+ = H2f{r) + f{r)Hi again in (47) to get instead 
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(El - El){n2je\f{r)\n,je) = {n2je\- - c'f"{r) - 2c'f'{r)^ ^^^^ 

- 2icarf'{r)V{r)\nije). 
The equations equivalent of Eqs. (31-33) are in this case 

{E, + E,){n2je\{-iarf{r))\n,je) = {n,je\- - c/'(r) - 2cf{r)^ ^^^^ 

- 2iarV{r)f{r)\nije), 

and 



cA+ 

{E2 - Ei){n2 j e I (-ia,/(r))|ni je) = {n2je\- cf{r) + -^(5f{r)+ ^^^^ 

2iarl3c^f{r)\nij€). 

We also have, for the matrix elements of (3f{r), 

{E2+Ei){n2je\(3f{r)\nije) = (^2 i e | - ica,/3/'(r)+ 

2[c^ + PV{r)] f{r)\mje). 



(52) 



These expressions are the basic equations for the case A^^ = 0. 

We can now obtain a recurrence relation valid in the diagonal case. First, let 
us use /(r) = in Eq. (48) to get 



{E2-Eif{n2je\r^\nije) =\{n2je\-c\X-iy-^ + c''^Pr^-^+ ^^^^ 

2ic^arpr^~^\nij e). 
Evaluating now equation (52) with /(r) = r^, we obtain 



{E2 + Ei){n2 j e \(3r^\ni j e ) ={n2 je\- icarl3Xr 



,2 ZP\ ^. (54) 



+ 2\^c'- -yj r^\nije), 
and eliminating the icar(3Xr^~^ between Eqs. (53) and (54), we finally get 



A+ 

[(^2 - Elf - Ac"] {n2 j e |r^|ni j e) = X^^{n2 j e |/?r^-2|ni 3 e ) 

- 4Zc2(n2 j e |/3r^- Vi 3 e) - 2c\E2 + Si)(n2 j e |/3r^|ni j e ) 
-c2A(A- I)(n2ie|r^-Viie). 
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This is the only recurrence relation we get in the diagonal case. To "close" the 
relation we can use the diagonal recurrence relations given in [9]. 
The special case when A = is of particular interest 

[(^2 - Elf - 4c^] 5n,n, = -^Zc'{n2je\^\mje) 

r (56) 

-2c\E2 + Ei){n2je\P\nije). 

This expression could be considered as a relativistic generalization of the Pasternak- 
Sternheimer rule of non relativistic quantum mechanics (Equation (8) of section II) 
[21], which says that the expectation value between hydrogenic states of the 1/r^ 
potential, vanishes when the orbital angular momenta of the states 1 and 2 coincide, 
i.e. when li — l2- In the relativistic case the expectation value of the /3/r potential 
(which could be regarded as the square root of 1/r^ including both signs), does not 
necessarily vanish even when the total angular momenta of the two states coincide: 
i.e. it does not vanish when ji = j2. Again, this agrees with the fact that the non- 
relativistic Paster nack-Sternheimer rule is applicable to eigenfunctions of potentials 
whose energy eigenvalues depend only on the principal quantum number — which is 
not the case for the hydrogen atom in Dirac relativistic quantum mechanics [14]. 
Moreover, two special cases are immediately deduced from this last expression 

(56): 

1) The first case, when ni ^ ^2, is 

Z(3 1 

{n2je\ |ni je) = --{E2 + Ei){n2j e\(3\nij e). (57) 

2) The other case follows when n\ = n2 

c^ = - {pv{r)) + E{P) = z(^^^+E{^), (58) 

which is the relativistic virial theorem [22]; from the relation < (3 >= E [9], we 
can also put it in the alternative form 



= ipVir)) +c' = -c'Z{^)+c'. (59) 



r 



V. The values of < r-^ > and < /?r^ >. 

The recurrence relations found above, involve in principle simple expressions 
(since they involve only matrix elements of Dirac hydrogenic states) that can be 
burdensome to handle. Given such situation, we have also calculated explicit formu- 
las that are needed to evaluate the diagonal and the non-diagonal matrix elements 
of interest. The expressions are related to the hypergeometric function and can be 
deduced from the two differential equations that follow from the Hamiltonian (20), 
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as it is shown in the Appendix. In particular, from Eq. (A. 15) we calculate < > 
and < (3r^ >. We quote the results here and refer to the Appendix for the details. 



< r"^ >= 



(60a) 



and 



< >= (2fc)A+iUi [^-(^)^' + In-in-iWv' + mc' «<;_i(A)] ; (606) 
in these expressions n = 0, 1, 2, • • •, and [23,27] 

(61) 

/ TflCp' — E I 

where ap — 1/137 is the fine structure constant and the /^^(A) symbols are defined 
in equation (A. 15) of the Appendix. The numbers u and v are constants such that 



u= {tj +s + n- Cz^"^)^^^, V = (n + 2s){Tj +s + n- Ci^"^)"^/^; 



(62) 



in the Appendix we give a simple proof of this result. Notice that in this section 
we have explicitly written h, e, and c in our results. Finally, to obtain C, we use 
relations (61) to get (tj + s + n — C^~^)~^ = (n + s — tj — C,v~^)/n{n + 2s); we 
need also (n + s) = C,E /yJmP'C^ — E'^, which is obtained from the expression for the 
energy eigenvalues of the Dirac hydrogen atom: 



-1/2 



E = mc^ 



\ (n-i - 1/2 + V(i + 1/2)2 -Z2a^)" 
elementary algebra gives then the result 



\C\ = 



n\ k 



[r(n + 2s + l)]-'/', 



Zac^ V 2m^ 

where we have written explicitly the dimensional factors. 



(63) 



(64) 
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Appendix. Explicit expressions for relativistic matrix elements of and 

It is possible to obtain explicit expressions for the diagonal and non diagonal 
matrix elements in the case V{r) = —Z/r. The purpose of this appendix is to give the 
basic relation that is needed for such evaluation. As we heavily draw from results 
previously obtained, in this section we use the notation of Ref. 23, in particular 
h = c = e = 1, though we sometimes write all the dimensional constants. 

We are interested in the bound states of the problem, so the quantity k = 
V m? — is positive. We can write the differential equations for the radial part of 
any central problem in terms of the dimensionless variable p = kr [23,25] and the 
symbols defined in (60) 

where we look for solutions of the form 

F{p) = V^^^ [V^_(p) + Mp)] , {A2) 



(Al) 



G{p) = x^T^ [i^- {p) - (p)] . (^3) 

The solution to these coupled differential equations can be written in terms of 
the Laguerre polynomials of non-integer index [25,27,28] 

V+(p) =ap^ exp(-p)£2^_i(2p), 
V^_(p)=6p^exp(-p)£^^(2p), 

where the Laguerre polynomials Ci^{p) are related to both the hypergeometric fun- 
ction, n, q; + 1; p), and the Sonine polynomials, Ta^\p) [28], through the rela- 
tion 
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^^^^^ = ^IrrlaVl^ ^^^"^' « + !;/>) = (-l)"r(a + n + l) ri-)(p), 

and a and b are constants. Substitution of these results in Eq. (Al) gives the condi- 
tion 



afr,- +s- Cf + n) + bin + 2s) = 0, 

(A6) 

6(rj — s + Ci' — n) — ari = 0. 

Solving these last two equations give us a relationship between n and u. Prom Eq. 
(45) we see that we can solve for the energy E and obtain the relativistic energy 
spectrum (63), provided we first introduce the principal quantum number N = 
j + 1/2 + n. To proceed further, we take 

}) = -a{Tj + s + n- Ciy~'^)/{n + 2s), {A7) 
and write the result in a symmetrized form: 



F{p) = Vmc^+ECp'e-P [u Cl\2p) + v ClUi'^p)] , 
G{p) = -Vmc^-ECp^e-P [uCl\2p) - v ClU{2p)] 

where 



(^8) 



C is a normalization constant that can be obtained from 

Jo ^• 

after some work we obtain 

i'^i = lJ/SF("+^'+^)r-'^ (-411) 

We can also calculate the expectation values for diagonal and non diagonal 
matrix elements. For diagonal, arbitrary power matrix elements of the form < > 
and < /3r-^ >, we need to calculate the expression 

/•OO 

InmW= e-^x'^+^C^ix)C^ix)dx. (A12) 
JO 

This expression converges for Re{a + A + 1) > 0, and is zero if A is an integer such 
that m — n > X> 0, where without loss of generality, we assume that m > n. Prom 
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Rodrigues formula and {dJ^ / dx'^)x^'^^ = (— — X]jnx'''^^ where [n], n an 
integer, is a Pochhammer symbol [28], we find, after a m-times partial integration, 



ra (^^- 1 \-( ,,fc r(n. + n + l)r(n + A- + A + !)[-/. -A]^ 
^nmW-—,2^^ 1) /,! (n-/c)!r(a + fc + l) " ^^^"^^ 

We use now the identity [—k — A]^ = [—k — A]fc[— A]^_fc, change the order of sum- 
mation k ^ n — k and use the identities 



[~X]m-n+k — [~X]m-n[~X + 171 — n]k, 

T{n + « + !) = (-1)'' T{a + n - k + 1) [-a - n]k, 
r(a + A + +n + 1) = (-l)''r(a + A + n- fc+l)[-a-A - n]k, (^^^) 

r(A+l) [-A-nJfc 

to obtain that 

r(a + A + n + l)r(A + n + l) 
^"'^^"^-'^ m!n! r(A + l) ^^^^^ 

3F2(— a — n, —A + m — n, —n; —A — n, —a — A — n; 1). 

We consider two cases for the general matrix elements (77.2^2 ^2\r^\niji ei) and 
{n2 32 ^2\(^r^\'n\ ji ei); the first one, when ki = k2, where we need to evaluate 

rOO 

K'n' (A) = / x^^+^^+^e--42^^) {x)£g^^^ (x) dr, (AW) 
Jo 

and the second one, when ki 7^ /c2, where we need 

K^l^^{X)= / r''+''+^e-^'''+''^^'"C^^''\2kir)Cg''\2k2r)dr. {A17) 
Jo 

In the first case, we see that integral (A16) is convergent if Re{si + S2 + A + 1) > 0, 
and vanishes when si — S2 + A is an integer such that m — n > si — S2 + A > 0. Using 
a similar reasoning as in the diagonal case, we get 

r , 1 r(si + S2 + +A + n + l)r(A + si-S2 + n+l) 

-A + S2 - S-i_\m-n , , T^TT 

ml n\ r(A + si — S2 + 1) 

3F2(-2si - n, -X + S2 - si + m - n, -n; -A + S2 - si - n, -X - Si - S2 - n; 1). 

(A18) 
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In the second case, the integral converges for Re{si + S2 + A + 1) > 0, and is 
not zero provided fci ^ ^2. A straightforward calculation by parts shows that 



nmK^) 2^2^ i\j\(m-i)\in- 3)\ 

r^,^^^.^,) ^(m + .2 + .l + A-^ + l), 

where ki ^ k2- Although less practical than the other expressions found here, we 
still can rewrite Eq. (A19) in a different form using the well known identities [28] 

bU-. = (-ir- ^^-"^'^ 



[—p — m + l]i T{—p — m + 1) ' 
r(— p — m + z + 1) = [—p — m + l]ir(—p — m + 1), 

mi{-iy 

{m-iy. = — 

[-m\i 

T{p +1) = [p- m+ l]rn^(,P - m+1). 

for any number p and m and i integers. After some algebra, we finally get 



(-1)"T(/; + 2.S1 + l)r(.si + S2 + A + 1) (-l)'^' [.si - ,s-2 + A + j - m + 1] 



E 



m!(fc2 + fci)*^+*^+^+^ ^ j\in- j)\T{2si+j + 1) 

k2 — k\ 

X 2Fii—m, si + S2 + A + 1; si - S2 + A + j - m + 1; — ). 

K2 + Kl 



{A21) 
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